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ABSTRACT 

Aims. We aim to better understand the heating of the gas by observing the prominent gas cooling line [Cii] at 158//m in the low- 
metallicity environment of the Local Group spiral galaxy M 33 at scales of 280 pc. In particular, we aim at describing the variation of 
the photoelectric heating efficiency with galactic environment. 

Methods. In this unbiased study, we used ISO/LWS [C ii] observations along the major axis of M 33, in combination with Herschel 
PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA H i data to study the variation of velocity integrated intensities. The 
ratio of [C ii] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated 
regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. 
Results. The heating efficiency stays constant at 0.8% in the inner 4.5 kpc radius of the galaxy where it starts to increase to reach 
values of ~ 3% in the outskirts at about 6 kpc radial distance. The rise of efficiency is explained in the framework of PDR models 
by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the 
significant fraction of H i emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) 
traced by Hi emission, the CNM contributes ~ 15% to the observed [Cii] emission in the inner 2 kpc radius of M33. The CNM 
contribution remains largely undetermined in the south, while positions between 2 and 7.3 kpc radial distance in the north of M 33 
show a contribution of ~ 40% ± 20%. 

Key words. Galaxies: ISM, ISM: photon-dominated regions, structure, evolution 
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1. Introduction 

In photon dominated regions (PDRs), FUV photons from stars 
dominate the chemistry and the energy balance in the interstellar 
gas. All the atomic and a large part of the molecular hydrogen 
of the ISM are located in PDRs which emit an important frac 
tion of far infrared (FIR) and milUmeter emission (Tielens & 
Hollenbach|T985l [Bakes & Tielensj 19941 [Hollenbach & Tielens 



19971. 



The [Cii] FIR fine structure line at 157.7jum is the most im- 
portant gas coolant. The [Oi] 63fim fine structure line starts 
to dominate in dense r and warmer regio ns when densities ex- 
ceed about 10"^ cm"^ ( Rolhg et al.||2006 i. The photoelectric ef- 
fect provides one of the dominant gas heating processes in 
PDRs. FUV photons eject electrons from dust grains or PAH 
molecules, heating the gas with their kinetic energy. Theoretical 
models have predicted efficiencies epE of up to a few percent 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



( Weingartner & Draine|200T I, consistent with observations. The 
ratio of emerging [C ii] intensity over the infrared continuum ra- 
diated by the dust has often been used as measure of this effi- 
ciency. Observations of clouds in the MiUcy Way show variations 
over more than 2 orders of magnitude, between 10"'* and 3 10"^ 



£l (e.g.|Vastel et al.|2001[pabart et al.|200T]|Mizutani et al.|2004[ 



[Jakob et al.|2007| l. A similar variation is found in observations of 
external galaxies (e.g. 'Malhotra et al.|200T] [Rubin et al.|2009" i. 
The scatter has been attributed to changes in the mean charge of 
small grains and PAHs (Okada et al. 2013| l. However, the change 
of epE in low-metallicity environments like those encountered in 
the Magellanic Clouds or M 33 is not yet well understood (e.g. 
Israel & Maloney "2011). Interestingly, the efficiency drops for 
local ULIRGs (Luhman et al.„2003 , Gracia-Carpio et al.[[20TT) 
as well as for some ULIRGs at high redshifts (e.g. [Stacey et al.[ 
[20T0l[Cox7tal.|2011| l. 



M33 is a nearby galaxy located at 840 kpc distance 
1991|l. Its overall metallicity is about half- 



( Freedman et al. 

solar ( Magrini et al.|2010|l, only sfightly larger than that of the 



Large Magellanic Cloud (LMC) ( [Hunter et al.|2007| l. M 33 is an 
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Fig. 1 : Herschel SPIRE 250^(111 map of M33 ( [Xilouris et al. 
2012| l. Units are Jy/(18"beam). Circles mark the position and 
beam size of the ISO/LWS [C ii] observations along the major 
axis of M33. The ellipse delineates the galacto-centric distance 
of 4.5kpc, dividing the observed inner (Nl and SI) and outer 
(N2 and S2) regions. The white square marks the BCLMP 302 
Hii region. Coordinates are R.A. and Dec. (Eq. J2000). 



Sc galaxy exhibiting a prominent, flocculent spiral structure to- 
gether with an underlying extended diffuse component as seen 
e.g. in the 250 /^m map of M 33 conducted in the framework of 



the Herschel open time key project HerM33es (Kramer et al 
2010 1 (Fig. 1). M33 has a moderate inclination of 56°, allow- 
ing studies of the interstellar medium with low depth along the 
line-of-sight. Its proximity allows high spatial resolution studies. 

While previous studies have discussed [C ii] emission only at 
few selected positions in M33, there has been no systematic 
study describing the spatial variation of [Cii] emission in the 
disk of the galaxy, jffigdonet al. (2003) used ISO/LWS to study 
[C ii] and other far-infrared (FIR) emission lines and the contin- 
uum in the nucleus and six Hn regions. They found a range of 
[Cii]/F IRlws values of between 0.2% and 0.7%. [Brauher et aT] 
( p0dH] l compiled ISO/LWS data of 227 galaxies, including 23 
positions in M33 with [Cii] data. Plotting the [Cii]/FIR ratio 
for all galaxies they find variations between 10 '^ and somewhat 
more than l%. |Mookerjea et al.| ( [201 l| l and |Braine et al.| ( |2012| i 
analyzed the first Herschel/PACS and HIFI spectroscopic data 
sets of the HerM33es project finding [Cii]/FIR ratios between 
0.01% and 2% in a 2' x 2' region centered on the BCLMP 302 
H II region and a ratio of 1 . 1 % at the position of the BCLMP 69 1 
Hii region lying at galacto-centric distances of 2.1 and 3.3 kpc, 
respectively, along the major axis. 

In this paper, we present archival ISO/LWS [C ii] data along 
the major axis of M 33 out to 8 kpc distance. We study the radial 
distributions and correlations between [C ii] and the FIR contin- 
uum, CO, H I, and Ha. We also study the radial distribution of 
the [C ii]/FIR ratio. We compare the observations in M33 with 
data of star forming regions in the Milky Way, with other ex- 
ternal galaxies, including low-metallicity objects. Local volume 
densities and FUV fields of the [C ii] and CO emitting gas are es- 
timated using the Kaufman et al.| ( |T999| K99) PDR model. The 
observed [C ii] emission is also compared with an estimate of the 
[Cii] emission emitted by atomic clouds. 



2. Observations and data analysis 

2.1. [Cii]l5Eijm (ISO/LWS) 

We list all ISCfjLWSF] [C ii] spectra observed along the major 
axis of M 33 at a position angle of 23° in Tablefl] The spectra 
were observed using the partial grating scan mode (LWS AOTjj 
L02). This AOT covers the wavelength range 43 - 196.9 fim and 
it has a medium spectral resolution of AA/A ~ 200 correspond- 
ing to 1500 km s ' at the wavelength of the [Cii] line. The LWS 
flux calibration and its relative spectral response function were 
derived from observations of Uranus (Swinyard et al.||1998]l. 
The angular resolution is 69.4" (Gry et al. 2003) which cor- 
responds to a linear resolution of 280 pc. The spectra had been 
automatically processed by the ISO system, passing scientific 
validation. We retrieved spectra at 77 positions from the ISO 
Data Archive (IDA) for further processing. The observed posi- 
tions cover a range of about + 8 kpc (±33') from the nucleus on 
a grid of about 208 pc (Figure [Til. 

We averaged the observations at each position, subtracted lin- 
ear baselines, and fitted a Gaussian to the line. Data were an- 
alyzed using the ISO Spectral Analysis Package (ISAP v2.1 



al.|1998[). In the follow ing, we assume a calibration er- 



[Sturm et i 

ror of 15% ("Higdon et al. 2003| l. Some sample spectra are shown 
in Appendix C] 

Figurel2]shows the variation of [C ii] intensities along the ma- 
jor axis. [Cii] is detected above 3cr at 36 positions. Stacking 
of neighboring positions increased the number of detections. 
Horizontal errorbars indicate the region over which the data were 



stacked. Tables A. 1 A. 2 in the Appendix list the stacked posi- 
tions. 



2.2. Far-infrared continuum 

To measure the total FIR continuum, we combined SPIRE and 
PACS maps of M33 at five wavelengths between 500 and 
lOOjum wavelength ( Boquien et al.|2011 Xilouris e t al.|2012) l, 
taken in the framework of HerM33es, with MlPS/Spitzer 24 and 
70;um maps ( Verley et al. 2007 ; Tabatabaei et al. 200T). 

These maps were smoothed to the ISO/LWS resolution us- 
ing Gaussian kernels. The fluxes were extracted using circular 
apertures of 69.4" centered at the ISO/LWS positions. A two- 
component grey body function was fitted to the spectral energy 
distribution (SED) at each position following the method de- 
scribed in Kramer et al. pOTOj). Integrating between 42.5 yum and 
122.5;um ( Dale & Helou|20()2[ l yields the FIR surface brightness. 



The dust emissivity index y6 was fixed at 1.5, which was found to 
be the best-fitting value for M 33 (|Kramer et al.|2010[ [Xilouris 
et al. 2012). The total FIR luminosity is, however, robust against 
changes of/?. Figurel2] shows the variation of relative FIR inten- 
sities along the major axis (cf. Table B.l i. A few sample SEDs 
are shown in Appendix[C| 

We also integrated the fitted SEDs over the range 3-1000 jum 



wavelengths to estimate total infrared (TIR) intensities ( Dale & 
|Helou|2002[ ), thereby deriving the ratios between TIR over FIR 
which lie between 2.7 and 1.3 (Table [bT| . 

2.3. CO, Hi, Ha 

Complementary CO and Hi data were used as tracers of the 
molecular and atomic gas to compare with [C ii] emission. 



Infrared Space Observatory ( Kessler et al. 



Long Wavelength Spectrometer in ISO (|Gry et al. 
Astronomical Observing Template 



19961 



20031 
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Table 1 : Observation Log of all [C ii] observations of ISO/LWS along the major axis of M 33 



IDA Name Abbrev." TDT' (Obs. #ID) 



R.A/ 



Dec' 



Observer ID Ref. Type AOT 



M33S2 


S2 


59901107 


01h33m08.5s 


+30dl7m00.0s 


KMOCHIZU 


- 


Raster 


L02 


M33S 


SI 


78600403 


01h33m37.1s 


+30d31m34.7s 


KMOCHIZU 


(1) 


Raster 


L02 


M 33 Nucleus 




80800367 


01h33m50.9s 


+30d39m36.8s 


HSMITH 


(1,2) 


Point 


L02 


M33N 


Nl 


78600801 


01h34m07.3s 


+30d46m55.7s 


KMOCHIZU 


- 


Raster 


L02 


M33N2 


N2 


59900605 


01h34m36.3s 


+31d01m29.2s 


KMOCHIZU 


- 


Raster 


L02 



'"' Abbreviation used in the text. '*' Target Dedicated Time. *''' The given coordinates are the central position of each raster strip of 19 positions. 
References. (l) |Brauher et"arip008l i. (2) |Higdon etlTlpOOS) . 



The CO 2-1 Une was mapped with the IRAM 30m telescope 
by Gardan et al. (2007 1 and Gratier et al. ( 20101 1. These maps 
cover the major axis out to a distance of 8.5 kpc in the north 
and 6.5 kpc in the south. The CO map covers all ISO/LWS [C n] 
positions but the eight southern-most positions. 

We determined 3 cr upper limits of the integrated intensities 
using cr - VA^Avres T^^ with the number of channels A^ over 
the velocity extent of the line, the velocity resolution Avies and 
corresponding baseline rms 7"™|5^ 

The Hi VLA map of M33 ( jGratier et al.|[20T0| covers the 
entire galaxy out to 8.5 kpc radial distance. Hi is detected at 
all ISO/LWS positions. While no single dish data were com- 
bined with the interferometric observations, the total flux recov- 
ered over the entire galaxy by the interferometric observations 
alone corresponds to more than 90% of the flux measured at the 
Arecibo single dish telescope (Putman et al.|2009 1. 



We also used a map of Ha emission presented in Hoopes & 
Walterbos|(2000| l and by [Verley et al[1 ( |2007| l. These data were 
obtained at the 0.6 m Burrell-Schmidt telescope at Kitt Peak 
National Observatory (KPNO). Ho- is detected at 42 ISO/LWS 
positions. We have calculated upper limits by measuring rms lev- 
els of the map in regions without Ha emission close to ISO/LWS 
[C ii] positions. 

Intensities have been calculated by smoothing all data to the 
angular resolution of the LWS [Cii] data as described in the 
Appendix[B] 

The variation of relative intensities of CO, H i, and Ha is 
shown in Figurel2]and absolute intensities are listed in Table B. 1 



We stacked CO, Hi, and Ha over the same positions as [Cn 
(cf.Tables [Al|[A2l ). 



3. Results 

3.1. Correlation between [Cn] and FIR, Ha, CO, Hi 

The emission of [Cii], the FIR continuum, and Ha, are all 
well correlated along the major axis, especially in the inner 10' 
(Fig. l2]l. These tracers of star formation all peak at the nucleus 
and drop by more than one order of magnitude beyond ~ 20' 
radial distance. Closer inspection shows that the FIR continuum 
drops steeper than [C ii], a finding which is discussed further be- 
low. In general however, [Cn], FIR, and Ha all trace the spiral 
arms as well as the inter-arm regions. The close correlation is 
even more clearly seen in Figurel3] where we plot the various 
tracers against [C ii] emission. A close correlation of [C ii] with 
other tracers of star formation is well known for other sources, 
and has been seen e.g. also in the ISO/LWS maps of a portion of 
the northern arm of M 3 1 by Rodriguez-Fernandez et al. ( 2006| l. 
The atomic and molecular gas traced by H i and CO also drops 
with radial distance, however the drop is much steeper for CO 
than for H i (Fig. |2]i. In the inner part of the galaxy within 10' 
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Fig. 2: [C ii], FIR, Ha, CO 2-1, and H i normalized integrated in- 
tensities along the major axis of M 33. Closed symbols show data 
within 4.5 kpc of the nucleus (Nl, SI) and open symbols show 
observations in the outer galaxy (N2, S2). Horizontal error-bars 
show the region over which different [C ii] spectra from neigh- 
boring positions were averaged. 



radial distance, the distribution of CO and H i is rather flat, in 
contrast to that of e.g. the FIR continuum. The distribution of 
emission of both gas tracers is not symmetrical with respect to 
the nucleus. Instead, CO and H i peak near 10' (~ 2.5 kpc) to the 
north and only show a secondary maximum at the nucleus. The 
absolute maximum corresponds to GMC91 (Engargiola et al. 



2003] [Gratier et al.|2010)|Buchbender et al.|26l3| l and to cloud 



245 in [Gratier etal. I ( |2012| ) 
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Fig. 3: Correlation of [Cii] intensities with those of the FIR con- Fig. 4: Normalized radial distributions of the [Cii]/CO 1-0, 
tinuum, Ha emission, CO, and Hi in double-logarithmic plots. FIR/CO 1-0, and Hi/CO 1-0 intensity ratios on the erg-scale. 
All intensities are given in units of ergs"' cm"^ sr"^ Straight 
lines delineate the results of linear least squares fits to the data. 
Closed and open symbols distinguish between positions in the 
inner (SI, Nl) and outer (S2, N2) disk of M 33, respectively. 

Table 2: Linear least squares fits to the correlations between 
[Cii] and FIR, Ha, CO, Hi (cf. Fig. ^ of the form logFIR= 
a + b log[C ii] with correlation coefficients r. 



FIR- [Cii] Hq'-[Cii] CO(2-1)-[Cii] Hi- [Cii] 



a 


3.71 


1.08 


-1.52 


-9.74 


b 


1.33 


1.23 


1.27 


041 


r 


0.90 


0.79 


0.59 


0.30 



Figure [3]shows correlations between [Cii] emission and FIR, 
Ha, CO, and Hi. Linear fits are weighted by the errors along 
both axes which we assume to be 15% for [Cii] and Hi, and 
20% for CO and the FIR continuum. The fits confirm that [C n] 
is strongly correlated with FIR and Ha, with linear correlation 
coefficients (r) 0.90 and 0.79, respectively (Table[2]). In contrast, 
the correlation with CO is much weaker, ?-=0.59, and very poor 
for Hi, r=0.30. This conforms that [Cii] is a good tracer of the 
star formation rate at the ISO/LWS beam size scale. Points from 
the northern most region N2 deviate slightly from the fit both in 
the [Cii]-FIR and the [Cn]-Ha plots, showing a worse correla- 
tion between [C ii] and star forming tracers in the northern, outer 
galaxy. 

3.2. Radial variation of intensity ratios 

The radial distributions of the ratios [C ii]/CO 1-0, FIR/CO 1-0 
and Hi/CO 1-0 are shown in Figure 4] For this, we estimated 
CO I-O intensities from the CO 2-1 ine (cf. Appendix |B]). All 
three ratios show a minimum in the inner parts of the galaxy 
and an increase towards the outer parts. This reflects the steep 
drop of CO intensities seen already in the radial distribution of 
intensities (Fig.|2]). The Hi/CO ratio shows a minimum near the 
nucleus and rises steadily towards the outskirts over about two 
orders of magnitude. While this behavior is about symmetrical 
in the North and South of the galaxy, this symmetry is broken 
for the [C ii]/CO and FIR/CO ratios. The latter two ratios show a 
minimum at about 10' to the North near GMC 91 where the CO 
emission peaks (Fig.|2|. To the south of this minimum, the two 




Fig. 5: Normalized radial distribution of the [Cii]/Hi intensity 
ratio. 



ratios steadily increase. Towards the north, the two ratios also 
increase, but with larger scatter 

In general, the [Cii] emission drops more steeply with 
galacto-centric radius ffian the H i emission. This is most clearly 
seen in the northern part of the strip (Fig.lSll. In contrast, the 
southern outer disk shows strongly varying [C ii]/H i ratios ow- 
ing to the strong variability of H i emission (cf. Fig.l2]i. 

3.3. [Cii]/FIR ratio 

In the inner part of the galaxy (N 1 , S 1 ), the [C ii]/FIR ratio stays 
constant within the measuring accuracy at about 0.8% + 0.2% 
(Fig.|6]), but rises significantly and steeply in the outskirts to val- 
ues of ~ 3%. The rise of the [C ii]/FIR ratio is caused by the steep 
drop of FIR emission relative to the [C ii] (Fig.|2|i, as is also seen 
in a plot of [C ii]/FIR versus FIR (Fig. [7]). The abrupt increase of 
the [C ii]/FIR ratio at about 4.5 kpc radial distance, both in the 
north and in the south of the major axis, occurs where the mor- 
phology of the whole galaxy seen in the optical changes (e.g. 
Sharma et al.||20Il ), i.e. just beyond the location of the promi- 
nent spiral arms. It also occurs in a region where the Hi/CO 
ratio rapidly increases (Fig.|4]i. As the O/H abundance gradient 
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Fig. 6; Radial distribution of the [C n]/FIR ratio along the major 
axis of M 33. 
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Fig. 7: [Cn]/FIR ratio as function of LpiR. The blue solid line 
shows a linear fit to the whole data set of M 33 while the green 
dotted line shows a linear fit to N2 and S2 points of the outer 
galaxy only. The black dashed-dotted horizontal line shows the 
lower [Cii]/FIR value found in our data set of M33. The blue 
dotted horizontal line shows the average [C n]/FIR value in nor- 
mal galaxies (Malhotra et al.|200 1 , MOl). We also show data of 
Milky Way regions (,Stacey et al. 1991, S91) and the low metal- 



licity objects LMC, SMC, IC 10 (Israel etal. 



& Maloney 201 1 111), (Madden et al. 1997 



1996 196), (Israel 



M97). 



in M 33 is shallow with a slope of only 0.035 dexkpc ' and 

with no signs of a break (Magrini et al. |2010 1, the metallicity 
cannot be key for the sudden increase of the [C ii]/FIR ratio. 

Though only few studies of the radial variation of the 
[Cii]/FIR ratio in galaxies exist, observations of M31, M51, 
NGC6946, and the Milky Way indicate that in general this 
ratio rises with galacto-centric distance. See the discussion in 
Rodriguez-Fernandez et al. ( 2006, 1 who present ISO/LWS maps 
of a portion of the northern arm of M 3 1 at 1 2 kpc radial distance 
showing a rather constant and high ratio of 2% while the nucleus 
shows only 0.6%. As M 31 is at the same distance as M33, the 
LWS observations of M 3 1 sample the same linear scale as the 
present observations of M33. Further below, we attempt to in- 



terpret the radial variation found along the major axis of M 33 
using PDR models. 

In M33, the [C ii]/FIR ratio drops with increasing FIR lumi- 
nosity. A power law fit to the entire data set results in [C ii]/FIR oc 
L^i^* (r=0.83). A fit to the data of the outer disk only results in 
a slightly steeper slope (Fig. IT]). In this figure, we also compare 
the M33 data with those of other galaxies. The M33 data are 
located in the region of low FIR luminosities Lfir=10^-10^ L© 



relative to those of the normal galaxies observed by Malhotra] 



etal. ( 2001 1 and high [C n]/FIR ratios of between 0.6% and 3% 



These data he in the same region as the data of the LMC, SMC 
and ICIO. 

The far-infrared luminosities of the low-metallicity systems 
are not directly comparable as their distances diff'er. However, 
while the LMC and SMC data were taken at 14 to 16pc linear 
resolutions, respectively ( [Israel et~ar 1996 Israel & Maloney 
2011 ), the [C ii] and FIR values were taken from maps where the 
fluxes were averaged over individual sources of between 35 and 
70 pc in size. The IC 10 observations have a resolution of 290 pc 
( |Maddenetal.|1997 !), similar- to that of the M 33 ISO/LWS data. 

The M 33 data in Figure |7]lie in the high limit of [C ii]/FIR ra- 
tios found in a sample of normal galaxies studied with ISO/LWS 
by |Malhotra et al. ( 2001, ). The latter appear on the right side of 
the plot with [C ii]/FIR = 0.01% - 0.7% and LpiR = 10^- lO" L© 
showing an average [C n]/FIR ratio of 0.32%. 

The [C ii]/FIR ratios in M 33 are 25 - 200 times higher than in 
Milky Way star forming regions for a similar LpiR range, while 
this factor decreases an order of magnitude (2.5 - 20) in Galactic 
GMCs. 



3.4. M33: A bridge between dwarf and normal galaxies ? 

After comparing [C ii]/FIR differences between M 33 and normal 
galaxies we also compared CO/FIR and [C ii]/CO ratios found in 
M 33 with those of other sources (Fig. [8]). 

In these plots, we find four mostly disjunct groups of ratios. 
There is a group formed by low metallicity objects and M 33 of 
the outer galaxy (N2, S2). The M 33 points of the inner galaxy 
(N 1 , S 1 ) form another group, which seems to connect low metal- 
licity objects with normal galaxies, and Milky Way GMCs (S91) 
which form the third group. Finally, PDRs in the Milky Way ex- 
posed to high FUV fields (MW SF regions, S91) form the fourth 
group. For normal galaxies, the [Cii]/CO 1-0 ratio increases 
with [Cii]/FIR. For M33 and other low metallicity objects the 
FIR/CO 1-0 ratios increase strongly with increasing [C ii]/CO 
1-0, a correlation which is also seen for normal galaxies though 
with increased scatter. The data points of the inner disk of M 33 
lie in-between the data of the normal galaxies at low [C ii]/CO 
ratios and the data of the other low-metallicity systems showing 
high [C n]/CO ratios. 

On average, the inner parts of M33 (Nl, SI) show lower 
[C ii]/CO and FIR/CO ratios than the other low-metallicity ob- 
jects. This suggests that CO is less photo dissociated in the in- 
ner disk of M33 than in the outer regions (S2, N2) and in the 
other low-metallicity objects. This is further discussed in the 
next paragraphs when comparing the observed ratios with PDR 
models. 



4. Discussion 

4.1. Diagnostic diagram of [Cn]/FIR vs. CO/FIR 

In Figure l9] we plot luminosities of [C ii] versus CO, normal- 
ized with FIR luminosities. To this diagnostic diagram we added 
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Fig. 9: [Cii] versus CO, normalized with the FIR continuum. Big green filled circles show ISO/LWS data of the inner SI, Nl regions 
of M33 while open circles show data of the outer S2, N2 regions. Small green circles show PACS observations of the BCLMP302 
Hii region in M33 (Mil). In addition we show data from Milky Way GMCs and star forming regions, and from other galaxies, 
comp i led fro m [Stacey et al.|(1991 , S91), Hailey-Dunsheath et al.| ( |2010l HDIO), |Gr acia-Carpio et al. (201 1 GC11) (GC 11), [Israel 



etal. (1996 196), Israel & Maloney (2011 111), Madden et al. ( 1997 M97), Malhotra et al. (2001 MO 1), and Luhman et al 



( |2003| L03). The lowest [C ii]/CO ratio observed with ISO/LWS m M 33 is 1000 (lower blue dashed Ime), while the highest ratio 
is 41200 (upper blue dashed line). Black solid and dotted lines indicate lines of constant density n and FUV field Go, respectively, 
from the standard K99 FDR model with Ay = 10 mag and solar metallicity Z -\. The dashed-dotted green line shows a K99 FDR 
model result for Ay = 1 mag, Z - \,n - 10^ cm"^. The knee of this curve is for Go = lO"^ (cf. Fig. 18 in K99). 



Table 3: Range of densities and FUV fields consistent with the 
standard K99 FDR model for sources shown in the diagnostic 
plot of [Cn]/FIR vs. CO/FIR (Fig.|9]l. 



Objects 



[C ii]/(CO 1-0) (i[cm-3] 



M33(N1-S1) 


1000 - 8000 


< 10* 


< 


10' 


M33(N2-S2) 


5000-4.110'' 


— 




— 


MW-SF 


1000 - 5000 


10^-10* 


10^'- 


> 10* 


MW-GMCs 


500 - 1000 


lO'-lO* 


500-^ 


-> 103 


Normal galaxies 


500-4000 


lO^-lO* 


10^ 


-103 


ULIRGs 


500 - 4000 


10*-10« 


103 


-10* 


Starburst nuclei 


-4100 


lO'-lO* 


500- 


-> 103 


Non-starburst nuclei 


500-4000 


10^-10' 


10^ 


-103 


LMC, SMC, IC 10 


8000 - 10^^ 


— 




— 



Notes. For the outer regions of M 33 and for the LMC, SMC, and IC 10, 
the standard model fails. 



observations of other galaxies and the corresponding lines of 
constant FUV-fields and local volume densities of standard K99 
FDR models with a cloud optical extinction of Ay = 10 mag for 
a solar metallicity. For completeness, we show in the Appendix 
[D][C ii] vs. CO, but normalized to the total-infrared (TIR) lumi- 
nosities. 



4.1.1. Observations 

M33. The ISO/LWS [Cii]/FIR ratios observed on scales of 
280 pc in M33 vary between ~ 0.8% in the inner parts of M33, 
rising to up to 3% in the outer regions, as already presented in the 
previous sections. The [Cn]/CO 1-0 ratios vary between 1000 
and 41200 while the CO 1-0/FIR ratios vary between 4 10"^ and 
8 10-^ 

On scales of 50 pc (12"), the [Cii]/FIR ratio varies between 
0.01% and 3% over the 2' x 2' map of the BCLMF 302 H ii re- 
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Fig. 8: Diagnostic plots of ratios of [Cii], CO and the FIR con- 
tinuum for M 33 and four other types of objects, low-metallicity 
systems, normal galaxies. Milky Way star forming regions, 
and GMCs. Upper panel: [Cii]/CO 1-0 vs. [Cii]/FIR. Bottom 
panel: FIR/CO 1-0 vs. [Cn]/CO 1-0. All panels are shown on 
a logarithmic scale. Symbols used have the same meaning than 
those used in Figure l7| 



gion |Mookerjea et al.| ( |20TT] Mli:0 The bulk of [C ii]/FIR ratios 
lie in the range ~0.7%-l%. The CO 1-0/FIR ratios he between 
4 10"^ and 8 10 *=. Interestingly, the PACS observations of the 
2' X 2' BCLMP302 region at 2.1 kpc radial distance cover the 
same range of [C ii]/FIR and CO 1-0/FIR ratios as found with 
ISO/LWS along the entire major axis of M 33 out to 8 kpc. 

Other galaxies. The [C ii]/CO ratios found in the inner parts of 
M33 lie in the same range of values found in the bulk of the 
normal galaxies and ULIRGs shown here. The outer regions of 
M 33 show higher values, similar to those found in other low- 
metallicity systems. 

The ISO/LWS [Cn]/FIR ratios of M33 are higher than in 
normal galaxies which only exhibit [Cii]/FIR ratios of up to 
0.4%. The low metallicity galaxies LMC, SMC, IC 10 show high 



'' To convert the TIR values used by Mil into FIR, we used the 
FIR/TIR ratios derived from the greybody fits (Sec. |2.2[ Tab. |B.l| (. Two 
ISO/LWS positions, N49 and N50, lie in the region mapped by Mil. 
We took the average FIR/TIR value for the conversion. 



[C ii]/FIR ratios, comparable to those found in M 33, as already 
discussed above. 

Local ULIRGs show CO/FIR ratios of less than 10 ^ while 
normal galaxies show higher ratios of up to 5 10"^. M33 shows 
slightly higher peak ratios of up to 8 10"^. 

4.1.2. PDR model results. 

Observations consistent with the standard model. Most of 
the [C n]/FIR and CO/FIR ratios observed in the inner disk of 
M 33 lie in the parameter space of local densities and FUV fields 
spanned by the standard PDR of K99 (solar metallicity Z - 1.0, 
Ay - 10 mag, Fig.[9]l. These ratios indicate densities of less than 
~ 10"^ cm"^ and FUV fields of less than Go ~ 10^^ in units of the 
local interstellar value. Some of the ratios observed in the inner 
galaxy and all ratios observed in the outer disk, are however not 
consistent with this standard model, as discussed further below. 

In the literature, the observed CO intensities have sometimes 
been multiplied by a factor of 2 for the comparison with the 
Kaufman PDR model. Hailey-Dunsheath et al.| ( [20T0l ) argue that 
CO 1-0 is optically thick stemming only from the frong-side of 
FUV illuminated clouds while the FIR emission is in general op- 
tically thin stemming from the front and the backside of clouds. 
Here, we do not apply any factor but rather use the observed val- 
ues argueing that the optical depth of the galactic CO emission 
on scales of 280 pc will be reduced due to the velocity disper- 
sion by turbulence and large-scale gradients, and that therefore 
an ad-hoc factor of 2 does not seem appropriate. 

In the framework of the standard PDR model, normal galaxies 
and especially the ULIRGs, tend to be consistent with higher 
densities of up to a few times 10^ cm"^^ than typically found in 
M 33. The latter are also consistent with higher FUV fields of up 
to ~ lO'*. 

Table [3] summarizes the range of values returned by the stan- 
dard K99 PDR model for the whole data sample shown in 
Figurel9] Normal galaxy and starburst nuclei points with extreme 
high [C ii]/CO values were ignored here. 



Observations inconsistent with the standard model. All 
ISO/LWS [C ii]/FIR and CO/FIR ratios of the outer disk of M 33 
are inconsistent with the standard PDR model for any density 
and FUV field. In general, [Cii]/FIR ratios above 2% cannot 
be reproduced by the standard model for CO/FIR ratios of less 
than 10"^. For lower CO/FIR ratios, the upper limit of [C n]/FIR 
which can be modeled, drops. For instance, at CO/FIR~ 10"^, 
the upper limit of [C ii]/FIR lies near 0.5%. 

The scatter of [Cii]/FIR and CO/FIR data seen in M33 on 
the large scales sampled by the ISO/LWS beam is similar to 
the scatter seen on small scales of 50 pc in the small sub-region 
BCLMP302. However, variations of metallicities are not ex- 
pected on such small scales. The observed scatter of [Cii]/FIR 
and CO/FIR ratio must therefore reflect the variation of other 
properties of the emitting gas and dust. 

As shown in Figuresl9]and 10 K99 PDR models of low opti- 
cal extinctions allow to reproduce the high [C ii]/FIR ratios ob- 
served predominantly in the outer disk of M 33. Assuming a total 
column density con^esponding to an Ay of only 1 mag and solar- 
metallicity, the ratios observed in the outer disk are consistent 
with n ~ lO-' cm""* and Go ~ 1 while the ratios observed in 
the inner disk are consistent with somewhat hi gher densities and 
FUV fields of « ~ 10"* cm"^ and Go ~ 10 (Fig!|l0|. The ratio of 
Go/n is about lO"-' cm^ 
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Subsolar metallicities do not need to be invoked to ex- 
plain these data sets, as the curves for models with subsolar- 
metallicity Z = 0. 1 are similar to the curves for Z = 1 .0. This 
also shows that [Cii] is relatively insensitive to changes in the 
metallicity. When [C ii] dominates gas-cooling, the models indi- 
cate that the gas temperature of the PDR surface layer emitting 
[Cii] adjusts in a way that the emergent [Cii] flux equals the 
FUV flux times the gas heating efficiency ( |Kaufman et al. 2006 1. 

As the observed metallicity gradient of M33 is shaUow 
(Magrini et al. 20101, we do not suspect that the average op- 
tical extinction of clouds changes abruptly from about 10 mag 
to about 1 mag between the inner and outer disk of M33. We 
suggest instead, that the overall low-metallicity environment of 
M 33 is composed of clouds of in general low optical extinctions 
of typically only 1 or a few magnitudes. 

In this scenario, the rather abrupt increase of heating effi- 
ciency beyond about 4.5 kpc radial distance is caused by a drop 
of average local densities of the molecular gas and a decrease of 
FUV field strengths. 

This drop of the FUV along the major axis of M 33 is al- 
ready indicated by the radially averaged extinction corrected 
FUV fluxes (cf. Fig. 3 in |Verley et al.|2009| l. [\^ey et al.| ( p007] l 
did IR photometry of 515 compact sources. The best models to 
reproduce the extinction seen in these H ii regions are the ones 
with Ay < 10 mag. 

The change of density from 10'* in the inner disk to lO"* in the 
outer disk, would imply that PDRs in the outer disk are typically 
about a factor 10 larger in size than in the inner disk: 
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Fig. 10: [Cii] versus CO, normalized with the FIR continuum 
(Figure 18 of K99). Black solid and dotted lines indicate lines of 
constant density for metallicities of Z = 1.0 and Z = 0.1, respec- 
tivey, from the K99 PDR model with Ay = 1 mag. Big green 
filled circles show ISO/LWS data of the inner SI, Nl regions of 
M 33 while open circles show data of the outer S2, N2 regions. 
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We estimated the FUV field in Habing units ( |Habing|ri968|l 
also directly from the FIR continuum, us i ng Gn=4;rlFfg/1.6-10~ 



in erg s cm "^ ( ,Mookerjea et al.||2011t [Kaufman et al.||1999| l. 



obtaining Go = 28 for the nucleus and Go = 1.3 for the N66 
position in the outer disk of M 33. The low FUV field found in 
the outer disk is consistent with the radiation field predicted by 
the K99 PDR model from flie [C ii]/FIR and CO/FIR ratios for 
low extinctions. 

For a fixed CO/FIR ratio, and a given density and metallic- 
ity, the modelled [C n]/FIR decreases with increasing optical ex- 
tinction (Figs. 16-18 in Kaufman et al. 1999 1. For example, at 
log(CO/FIR)=6, log([Cii]/FIR) decreases by a factor ~ 5 from 
-1.5 at 1 mag to -2.2 at lOmag. In clouds of low optical extinc- 
tions, FUV photons penetrate deeper into the cloud. For spheri- 
cal clouds, the smaller CO core is surrounded by a larger [C n] 
emitting region, leading to enhanced [C ii]/CO ratios (cf. IBolatto 
[etTl . 1999 ; R olhg et aipOOel . The high observed [C ii]/FIR and 
[C ii]/CO ratios are best explained by clouds of low columns (op- 
tical extinctions) leading to an increase of the [C ii] layer rela- 
tive to the total gas traced by CO, and also relative to the to- 
tal dust traced by the FIR continuum. We propose that reduced 
volume densities and the geometrical dilution of the FUV field 
when photons can enter more deeply, are of secondary impor- 
tance only. 

The low optical extinctions are consistent with a reduced dust- 
to-gas ratio in the low metallicity environment of M 33. As dis- 
cussed by Israel et al. ( 1996 ) and Israel & Maloney (201 1 ) using 
data of the LMC, the ratio of total gas column density over opti- 
cal extinction, Nh/Av, increases in low metallicity environments. 




Dista.nce [kpc] 



Fig. 11: Fraction of Hi column density stemming from PDRs. 
EiTorbars only show the 15% observational error of Hi in- 
tensities. Two positions in the south for which A^(Hi,PDR)> 
A^(H i,obs) are not shown. 



4.2. Hi emission from PDRs 

In the following, we will assume that all observed Hi emis- 
sion stems from the cold neutral and atomic medium (CNM) 
and from photon-dominated regions. Before discussing a possi- 
ble contribution of [C ii] emission from the CNM, we first need 
to study the Hi fraction stemming from PDRs. At the surface of 
photon dominated regions, molecular hydrogen is photodissoci- 
ated and a layer of atomic hydrogen is formed. The column from 
H I is a function of impinging FUV field over the volume density 
of the molecular gas, Go/n, as shown e.g. by Sternberg ( 1988|. 
|Stacey etal.| ( |199l) studied the contribution of H i emission from 
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PDRs to the total hydrogen column in a sample of galaxies. With 
the typical densities and FUV fields derived above for M 33, we 
find a typical Go/« ratio of 10"^ cm^ in the inner and outer disk 
of M 33. Next we use Eq. 3 of |Heiner et aL] ( [20n] ): 



7.8 10-" , 
N(ni, PDR) = ln( 



(l-^^-'")' 



(2) 



with a dust-to-gas ratio D normalized to the solar neigh- 
borhou d value of (12 + log(0/H)) - 8.69 (Eq. 2 of [Heiner et al 
( 2011| l) and a constant oxygen abundance of 12 -H log(0/H) = 
8.27 (see references in |Buchbender et al.|2013| l. 

This leads to a constant Hi column density stemming from 
PDRs of 3.25 lO^^cm-^. 



In Figure 1 1 we compare the H i column density from PDRs 
with the observed beam-averaged H i column density, assuming 
that the PDRs fill the beam. We did not attempt detailed PDR 
modeling to derive the beam filling factor, e.g. by comparing the 
observed FIR continuum with the best fitting FUV field from 
models. (Stacey et W. ( 1991| l derive a typical beam filling factor 
of 0.3 for the sample of galaxies they studied. In M33, the es- 
timated fraction of H i column densities stemming from PDRs 
stays between 15% and 70% for all data between 2kpc radial 
distance in the south and 7.3 kpc in the north. Two positions fur- 
ther south exhibit very low H i columns which are lower than the 
estimated H i column from PDRs. Their A^(H i,PDR)/A^(H i,obs) 
ratios rise above 1 to values of up to 4, indicating that for these 
regions the underlying assumptions of Go, n, and the beam filling 
factor are not valid. 
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Fig. 12: Fraction of the observed [Cn] emission stemming from 
the cold neutral medium for T-SQ K and «=100 cm"^, after cor- 
rection of the observed Hi for the contribution from PDRs. 
Squares mark the individual data along the major axis. Their er- 
rorbars only include the 15% observational error of the H i inten- 
sities. Crosses in red mark the fractions averaged over the inner 
2 kpc, the northern 2 to 4.5 kpc, and the outer, northern disk be- 
tween 4.5 and 7.4 kpc. Corresponding errorbars show the stan- 
dard deviation of the individual data and the radius interval over 
which the values were binned. 



4.3. [Cn] from the atomic medium 

Here, we investigate which fraction of the observed [C ii] emis- 
sion could stem from the atomic cold neutral medium (CNM) in 
M33. 



As described in Section pTT] the emission of H i drops only 
weakly towards the outskirts, while the CO emission drops 
steeply (Fig.l2]i. In the southern part of the strip, the Hi/CO rises 
continuously from a minimum in the center towards the outskirts 
at about 5 kpc radial distance (Fig. HI). In the North, on the other 
hand, the H i/CO shows an abrupt increase at about 4.5 kpc ra- 
dial distance, i.e. at the same distance where the [C ii]/FIR ratio 
rises to values of ~ 3% (Fig.l6]l. The relative contribution of the 
atomic gas tracing the CNM to the total gas emission (CO-i-H i) 
rises with galacto-centric distance. Also, the Hi/[Cii] intensity 
ratio rises strongly in the northern part of the strip over almost 
an order of magnitude (Fig.lSll. This is because the [Cii] emis- 
sion drops more steeply with radius than the H i. In the north, 
the observed [C ii] emission has dropped by a factor 27 relative 
to its value at the nucleus while Hi has dropped only by 40%. 
Any contribution of the CNM to the [C ii] emission relative to 
the contribution from PDRs should be most prominent in the 
outskirts. 

We note that the correlation between H i and [C ii] emission 
is rather weak throughout the strip, and also in the outskirts 
(cf. Fig.[3]l. pookerjea et al.| ( |20TT] l and |Braine elaTl l pOTIi) stud- 
ied [C ii] HIFI spectra in two regions of ~ 1' - 2' along the ma- 
jor axis of M33 at radial distances of 2.1 and 3.3 kpc. In both 
sources, the Hi lines are shifted systematically relative to the 
[C ii] lines by about 5 kms ' . In addition, H i linewidths are about 
30% broader than those of [Cii]. These findings indicate that 
[C ii] emission, at least in these two regions of the inner galaxy, 
does not simply trace only the atomic medium. 

To quantify a possible contribution of the cold neutral medium 
(CNM) to the [C ii] emission, we use the H i line intensity, cor- 
rected for the contribution from PDRs, to derive the intensity 
ofthe[Cn]line(cf . Crawf ord et al.|1985[ [Madden et al.|1993l 



1997 Langer et al.j2010). Here, we assume optically thin [Cii] 



and H i emission and beamfilling factors of unity. We estimate 
the [C ii] intensities stemming from the remaining fraction of the 
atomic gas via: 



/'(HI) = 3.35 10"'* /(HI) 
A^(HI) = 1.82 10"* /'(HI) 
A?(HI,CNM) = A^(HI) - A^(HI, PDR) 
A?(C+,CNM) = Xc. A^(HI, CNM) 
/(C+,CNM) = 2.35 10"2'A^(C+) X 

2 exp(-A£'/r) 



(r 



2 exp(-A£'/r) + («„/«) 



'-2 1 

+ (n„ln)l 



(3) 
(4) 
(5) 
(6) 
(7) 

(8) 



with intensities / in erg s"' cm"^ sr"', and /' in Kkms"' 
column densities in cm"^, the energy of the [C ii] ^P 
above ground AE - hv/k = 91.3 K, and the critical density 
for collisions with hydrogen atoms and/or molecules of «rr = 



, and 
3/2 State 



3 10 cm (cf. discussion and references in |Langer et al.|2010| l. 



The fractional abundance X(C^) of C^ in the gas with respect to 
H I li es in the range 1.4 to 1 .8 10 "* in the local ISM of the Milky 
Way ( [Sofia et al.|1997| l. For the average low metallicity enyiron- 
ment of M33, we assume X(C+)=0.6 10"'* (Magrini et al.|2010 



Henry et al. 2000 1. In addition, we assume diffuse atomic clouds 
of n = 100 cm~^ and T = 80 K corresponding to a pressure of 
about 8000 K cm"-' typical for the difiiise atomic medium in the 
Milky Way (e.g. [Wolfire et al.| 1995 [[Dickey et al.|2000| . Higher 
temperatures and densities would increase the [Cii] intensities 
but quickly become inconsistent with the generally accepted val- 
ues for the pressure of the Galactic atomic ISM (cf. discussion 
in Sec. 4.2.1. of [Madden et al.|1993| ). 
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Figure[T2] shows the estimated fraction of the observed [C n] 
intensity which stems from the CNM. As expected from the ob- 
served ratio of H i over [C ii] intensities, the fraction is lowest in 
the nucleus and in general rises towards the outskirts. Towards 
the southern half, the CNM fraction shows strong variability re- 
flecting the variations of Hi emission on small scales (cf. Fig.|2]l. 
However, towards the north, the fraction rises in general from 
20% ± 6% in the inner 2kpc, to 37% + 14% between 2 and 
4.5 kpc, to 43% + 18% between 4.5 and 7.4 kpc radial distance. 

The rise of [C ii]/FIR by a factor of more than 3 with distance 
may therefore be partially explained by a rising contribution of 
the CNM to the [C ii] emission while the [C ii]/FIR ratio of PDRs 
rises only moderately with only somewhat lowered values of Ay, 
n. Go relative to the PDR models for the inner galaxy of Ay ~ 
10 mag, n ~ 10"^ cm"^. Go ~ 10. A better knowledge of the 
temperature and density of the CNM is needed to derive more 
quantitative conclusions. 

Here, we conclude that PDRs dominate in the inner parts of 
M33 while both, the CNM and PDRs, can explain the [Cn] 
emission observed in the outskirts. 



4.4. The importance of the [Oi] 63jum line 

A more complete study of the major gas cooling lines of the ISM 
would need to include the [O i] 63yum line as a tracer of high den- 
sities in PDRs. Higdon et al. ( 2003| l detected this line at scales 
of 280 pc in the nucleus and towards five Hii regions in M33, 
finding [Cii]/[Oi] ratios between 1.6 and 2.2. Towards the Hii 
region BCLMP 302 and on scales of 50 pc, the ratio varies be- 
twen 2.5 and 10 ( |Mookerjeaetal.|2011| l. Further Herschel/PACS 
and HIFI observations of several regions along the major axis of 
M33 have been conducted in the framework of the HerM33es 
project, to obtain a more complete dataset which will allow e.g. 
to derive more accurate estimates of the photoelectric heating ef- 
ficiency and its variation on small scales, as also it will allow to 
determine gas properties (Ay, n, T) using PDR models. 



5. Summary and conclusions 

In this paper we have studied ISO/LWS [Cii] emission along 
the major axis of M 33 at scales of 280 pc. Stacking of few ad- 
jacent positions allowed to detect [Cii] between galacto-centric 
distances of ~ 6 kpc in the South to ~ 6.7 kpc in the North. 

The radial distributions of the [Cii]/CO and Hi/CO ratios 
show an increase from the inner regions to the outer regions. 
While the Hi/CO ratio has a minimum near the nucleus, the 
[Cii]/CO ratio shows a minimum at about 10' to the North, cor- 
responding to the global maximum of CO emission along the 
major axis. 

The radial distribution of [Cii] shows a strong correlation 
with star formation tracers, the FIR continuum and Ha emis- 
sion, closely following the spiral arm structure. The correlation 
becomes weaker in the outskirts of M33 beyond 4.5 kpc radial 
distance to the nucleus. 

The [Cn]/FIR ratio remains constant at ~ 0.8% in the cen- 
tral 4.5 kpc radius and then increases rapidly to values of up 
to 3%. The [Cii]/FIR ratio is a strong function of the FIR lu- 
minosity. The highest [C ii]/FIR ratios are found for the lowest 
FIR luminosities. Other low metallicity systems (LMC, SMC, 
IC 10) show similar [Cii]/FIR ratios as M33 which are higher 
than those found in the Milky Way and in normal galaxies. 

The variation of [C n]/FIR and CO/FIR ratios found in a small 



2' X 2' H II region of M 33 on small scales of 50 pc (Mookerjea 



|et al.| 2011| l is very similar to the variation found with the 
ISO/LWS observations. This is another indication that metallic- 
ity variations play a minor role in determining the [C ii]/FIR ra- 
tio. 

Diagnostic plots of [Cii]/FIR vs. CO/FIR or FIR/CO vs. 
[C ii]/CO show smooth transitions from normal galaxies to the 
positions in the inner 4.5 kpc radius of M33, to the outer 
parts of M 33 together with the other low-metallicity systems. 
Observations of the inner disk of M 33 serve as bridge between 
normal galaxies exhibiting low [C n]/FIR and low [C n]/CO ra- 
tios to the low-metallicity systems showing high [C ii]/FIR ratios 
and also high [C ii]/CO ratios. 

The relative lack of dust shielding in the low metallicity sys- 
tems enhances CO photo-dissociation, and, hence, increases the 
[Cii]/CO ratios. The standard K99 PDR models of 10 mag op- 
tical extinction fail to reproduce the high end of the observed 
[C ii]/FIR and [C ii]/CO ratios. However, models of low optical 
extinction expected in the low-metallicity environment of M 33 
do reproduce the observed ratios for a constant metallicity. In the 
framework of models of Ay = 1 mag, the variation of observed 
[C ii]/FIR ratios is caused by variations of the local densities 
dropping from about 10"* cm"^ in the inner disk to 10^ cm"^ in 
the outer parts. In addition, FUV field strengths reach ~ 10 times 
the average interstellar radiation field in the inner disk while 
dropping to only a few in the outskirts. High FUV fields tend 
to lower the photoelectric heating efficiency as grains and PAHs 
become positively charged (Okada et al. 2013). In contrast, the 
high [Cn]/FIR ratios observed in low metallicity galaxies can 
be caused by geometric dilution of low FUV fields in clouds of 
normal densities where Carbon stays ionized. 

In addition to PDRs, the atomic cold neutral medium (CNM) 
traced by Hi may contribute to the observed [Cii] emission. 
Therefore, the above scenario may need to be revised. We con- 
clude that the CNM, corrected for the H i emission from PDRs, 
contributes ~ 15% to the [Cii] observed in the inner radius of 
2 kpc. However, the CNM may contribute about 40% of the ob- 
served [Cii] in the outer, northern disk where Hi emission is 
much stronger relative to [C ii] . 

Acknowledgements. We thank S.Hailey-Dunsheath and J. Gracia-Carpio for 
providing us data shown in Figure|9] We thank an anonymous referee for in- 
sightful comments. 
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from the CO 2-1 data, we used a linear function of the 2-1/CO 
1-0 ratio dropping from 0.8 in the nucleus of M 33 to 0.5 at a 
galactocentric distance of 8.5 kpc ( Gratier et al.|2010 ). Next, we 
smoothed the map to the ISO/LWS resolution and converted the 
data to intensities in erg cm^^s^'sr"^ 

Both the H i data and the Ha data were also smoothed to the 
LWS resolution. The original resolutions of the H i data was 11" 
while we assumed an original pencil beam for the Ha data. 

Table B.l lists for all positions and stacked areas the resulting 
integrated intensities of [Cii], Hi, Ha, CO 2-1, 1-0, the FIR 
continuum, and the TIR/FIR ratio. 

Appendix C: Sample [Cii] ISO/LWS spectra and 
SEDs 

Here, we show four examples of the [C ii] spectra and SEDs ob- 
tained. Figure C.la shows the only spectrum with a baseline fit 
of order 1. Figure C.lc shows the spectrum from the nucleus, 
where the lowest rms and highest flux are reached. The highest 
rms and also the lowest flux peak of the four spectra is shown 



in position N61 (fig. C.lg i which belongs to the northern, outer 
N2 region. The SEDs show the strongest warm dust component 
in the nucleus (Fig. |C.ld[ ) while the weakest warm component 
is seen in N61 (Fig.|C.lh|. Cold dust components are similar in 



the inter-arm region (Fig. C.lbi, the nucleus, and in BCLM302 



(Fig. C.lf I. The outer points also show weaker cold dust emis- 



sion than the other regions. 



Appendix D: [Cii]/TIR vs. CO/TIR 



For completeness, we show in Figure [DT| the diagnostic plot of 
[C ii] vs. CO luminosities, normalized with the TIR luminosities 



derived from the two-greybody fits described in Section 2.2 



Appendix A: List of positions 



The Tables |A.1| and |A.2| list the 77 positions observed with 
ISO/LWS, their galacto-centric distance, and whether [C ii] was 
detected. 



Appendix B: Intensities at the observed ISO/LWS 
positions 

The IRAM 30m CO 2-1 data on the T"^ scale were converted into 
Ttnb temperatures using the forward efficiency Feff - 0.90 and 
the main beam efficiency Betr - 0.49: Tf„[, = (/^etr/^eff) T'^. The 
CO 2-1 map was first smoothed from the original 1 1" resolution 
to 22" using a Gaussian kernel. To derive CO 1-0 intensities 
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Table A. 1 : ISO/LWS positions observed in [C ii] along the southern part of the major axis of M 33. 



PosID 


RA 


DEC 


Ang. Dis ['] 


Lin. Dis [kpc] 


[C ii] detection ? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


SI 


01 


32:57 


+30:10:06 


-32.70 


-7.99 


S1-S8 


S2 


01 


32:58 


+30:10:52 


-31.80 


-7.77 


S1-S8 


S3 


01 


32:59 


+30:11:38 


-31.00 


-7.57 


S1-S8 


S4 


01 


33:01 


+30:12:24 


-30.10 


-7.35 


S1-S8 


S5 


01 


33:02 


+30:13:10 


-29.30 


-7.16 


S1-S8 


S6 


01 


33:03 


+30:13:56 


-28.40 


-6.94 


S1-S8 


S7 


01 


33:05 


+30:14:42 


-27.50 


-6.72 


S1-S8 


S8 


01 


33:06 


+30:15:28 


-26.70 


-6.52 


S1-S8 


S9 


01 


33:07 


+30:16:14 


-25.80 


-6.30 


S9-S14 


SIO 


01 


33:08 


+30:17:00 


-25.00 


-6.11 


S9-S14 


Sll 


01 


33:10 


+30:17:46 


-24.10 


-5.89 


S9-S14 


S12 


01 


33:11 


+30:18:32 


-23.30 


-5.69 


S9-S14 


S13 


01 


33:12 


+30:19:18 


-22.40 


-5.47 


S9-S14 


S14 


01 


33:14 


+30:20:04 


-21.60 


-5.28 


S9-S14 


S15 


01 


33:15 


+30:20:50 


-20.70 


-5.06 


S15-S19 


S16 


01 


33:16 


+30:21:36 


-19.90 


-4.86 


S15-S19 


S17 


01 


33:18 


+30:22:22 


-19.00 


-4.64 


S15-S19 


SIS 


01 


33:19 


+30:23:08 


-18.10 


-4.42 


S15-S19 


S19 


01 


33:20 


+30:23:54 


-17.30 


-4.23 


S15-S19 


S20 


01 


33:24 


+30:24:40 


-16.40 


-4.01 


S20-S21 


S21 


01 


33:25 


+30:25:27 


-15.60 


-3.81 


S20-S21 


S22 


01 


33:27 


+30:26:13 


-14.70 


-3.59 


S22-S23 


S23 


01 


33:28 


+30:26:59 


-13.90 


-3.40 


S22-S23 


S24 


01 


33:30 


+30:27:45 


-13.00 


-3.18 


S24-S26 


S25 


01 


33:31 


+30:28:31 


-12.20 


-2.98 


S24-S26 


S26 


01 


33:33 


+30:29:17 


-11.30 


-2.76 


S24-S26 


S27 


01 


33:34 


+30:30:03 


-10.50 


-2.57 


X 


S28 


01 


33:36 


+30:30:49 


-9.60 


-2.35 


X 


S29 


01 


33:37 


+30:31:35 


-8.75 


-2.14 


X 


S30 


01 


33:39 


+30:32:21 


-7.89 


-1.93 


X 


S31 


01 


33:40 


+30:33:07 


-7.04 


-1.72 


X 


S32 


01 


33:42 


+30:33:53 


-6.19 


-1.51 


X 


S33 


01 


33:43 


+30:34:39 


-5.33 


-1.30 


X 


S34 


01 


33:45 


+30:35:25 


-4.48 


-1.09 


X 


S35 


01 


33:46 


+30:36:11 


-3.63 


-0.89 


X 


S36 


01 


33:48 


+30:36:57 


-2.79 


-0.68 


X 


S37 


01 


33:49 


+30:37:43 


-1.95 


-0.48 


X 


S38 


01 


33:51 


+30:38:29 


-1.14 


-0.28 


X 



Notes. Columns (4) and (5) list galacto-centric distances. Negative values indicate southern positions. Column (6) indicates whether [C ii] emission 
was detected above 3(T at a given position (marked by "X") or were averaged over a given range of positions. 
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Table A.2: ISO/LWS positions observed in [C ii] along the northern part of the major axis of M 33. 



PosID 


RA (J2000) 


DEC (J2000) 


Ang. Dis ['] 


Lin. Dis [kpc] 


[C II] Detection? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


39 (Nucleus) 


01:33:51 


+30:39:37 


0.00 


0.00 


X 


N40 


01:33:54 


+30:40:01 


.80 


0.20 


X 


N41 


01:33:55 


+30:40:47 


1.58 


0.39 


X 


N42 


01:33:57 


+30:41:33 


2.42 


0.59 


X 


N43 


01:33:58 


+30:42:19 


3.26 


0.80 


X 


N44 


01:34:00 


+30:43:05 


4.11 


1.00 


X 


N45 


01:34:01 


+30:43:51 


4.96 


1.21 


X 


N46 


01:34:03 


+30:44:37 


5.81 


1.42 


X 


N47 


01:34:04 


+30:45:23 


6.67 


1.63 


X 


N48 


01:34:06 


+30:46:09 


7.52 


1.84 


X 


N49 


01:34:07 


+30:46:55 


8.38 


2.05 


X 


N50 


01:34:09 


+30:47:41 


9.23 


2.26 


X 


N51 


01:34:10 


+30:48:27 


10.1 


2.47 


X 


N52 


01:34:12 


+30:49:13 


10.9 


2.66 


X 


N53 


01:34:13 


+30:49:59 


11.8 


2.88 


X 


N54 


01:34:15 


+30:50:45 


12.6 


3.08 


X 


N55 


01:34:16 


+30:51:31 


13.5 


3.30 


X 


N56 


01:34:18 


+30:52:17 


14.4 


3.52 


X 


N57 


01:34:19 


+30:53:03 


15.2 


3.71 


N57-N60 


N58 


01:34:21 


+30:53:49 


16.1 


3.93 


N57-N60 


N59 


01:34:23 


+30:54:35 


16.9 


4.13 


N57-N60 


N60 


01:34:24 


+30:55:21 


17.8 


4.35 


N57-N60 


N61 


01:34:26 


+30:56:07 


18.7 


4.57 


X 


N62 


01:34:27 


+30:56:53 


19.5 


4.76 


X 


N63 


01:34:29 


+30:57:39 


20.4 


4.98 


X 


N64 


01:34:30 


+30:58:25 


21.2 


5.18 


- 


N65 


01:34:32 


+30:59:11 


22.1 


5.40 


X 


N66 


01:34:33 


+30:59:57 


22.9 


5.60 


X 


N67 


01:34:35 


+31:00:43 


23.8 


5.82 


N67-N68 


N68 


01:34:36 


+31:01:29 


24.6 


6.01 


N67-N68 


N69 


01:34:38 


+31:02:15 


25.5 


6.23 


X 


N70 


01:34:39 


+31:03:01 


26.4 


6.45 


N70-N72 


N71 


01:34:41 


+31:03:47 


27.2 


6.65 


N70-N72 


N72 


01:34:42 


+31:04:33 


28.1 


6.87 


N70-N72 


N73 


01:34:44 


+31:05:19 


28.9 


7.06 


N73-N77 


N74 


01:34:45 


+31:06:05 


29.8 


7.28 


N73-N77 


N75 


01:34:47 


+31:06:51 


30.6 


7.48 


N73-N77 


N76 


01:34:48 


+31:07:38 


31.5 


7.70 


N73-N77 


N77 


01:34:50 


+31:08:24 


32.3 


7.89 


N73-N77 



Notes. Columns (4) and (5) list galacto-centric distances. Column (6) indicates whether [C ii] emission was detected above 3cr at a given position 
(marked by "X"), was averaged over a given range of positions, or whether no emission was detected (marked by "-"). 
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Table B.l: Integrated intensities at the positions with [C ii] detections. 



Pos. ID 


[Cii] 


Hi 


Ha 


CO(2-l) 


CO 1-0 


FIR 


TIR/FIR 


S1-S8 


2.20e-06 


6.88e-13 


3.05e-07 


— 


— 


9.12e-05 


2.52e+00 


S9-S14 


1.28e-06 


1.37e-13 


3.11e-07 


1.18e-09 


8.32e-ll 


4.25e-05 


2.50e+00 


S15-S19 


1.26e-06 


1.73e-12 


8.68e-07 


3.13e-09 


2.39e-10 


4.17e-05 


2.47e+00 


S20-S21 


4.70e-06 


1.48e-12 


3.75e-06 


6.35e-09 


5.15e-10 


8.39e-04 


1.54e+00 


S22-S23 


3.75e-06 


1.86e-12 


1.20e-06 


7.06e-09 


5.87e-10 


4.61e-04 


1.89e+00 


S24-S26 


2.69e-06 


2.00e-13 


9.64e-07 


1.70e-09 


1.44e-10 


2.62e-04 


2.32e+00 


S27 


5.04e-06 


7.98e-13 


3.63e-06 


3.10e-09 


5.51e-10 


4.82e-04 


1.44e+00 


S28 


8.38e-06 


1.08e-12 


7.35e-06 


5.92e-09 


1.04e-09 


8.19e-04 


1.48e+00 


S29 


1.42e-05 


2.11e-12 


1.34e-05 


1.69e-08 


2.94e-09 


1.73e-03 


1.49e+00 


S30 


1.57e-05 


2.02e-12 


1.61e-05 


1.81e-08 


3.11e-09 


1.72e-03 


1.73e+00 


S31 


1.45e-05 


1.68e-12 


1.24e-05 


1.28e-08 


2.17e-09 


1.62e-03 


1.56e+00 


S32 


6.97e-06 


1.35e-12 


7.45e-06 


9.68e-09 


1.63e-09 


1.34e-03 


1.38e+00 


S33 


9.54e-06 


1.25e-12 


5.81e-06 


8.88e-09 


1.48e-09 


1.16e-03 


1.44e+00 


S34 


1.43e-05 


1.74e-12 


9.92e-06 


1.41e-08 


2.32e-09 


1.78e-03 


1.42e+00 


S35 


1.78e-05 


2.08e-12 


1.66e-05 


2.75e-08 


4.30e-09 


2.00e-03 


1.32e+00 


S36 


1.64e-05 


1.55e-12 


1.41e-05 


2.02e-08 


3.15e-09 


2.01e-03 


1.39e+00 


S37 


1.41e-05 


1.02e-12 


1.38e-05 


1.50e-08 


2.34e-09 


1.85e-03 


1.37e+00 


S38 


1.75e-05 


1.04e-12 


1.90e-05 


1.75e-08 


2.74e-09 


2.25e-03 


1.34e+00 


39 (Nucleus) 


3.27e-05 


1.92e-12 


2.77e-05 


3.41e-08 


5.34e-09 


3.62e-03 


1.34e+00 


N40 


2.36e-05 


2.01e-12 


1.52e-05 


2.85e-08 


4.45e-09 


2.78e-03 


1.33e+00 


N41 


1.45e-05 


2.28e-12 


1.03e-05 


3.17e-08 


4.96e-09 


2.45e-03 


1.39e+00 


N42 


1.34e-05 


2.09e-12 


8.94e-06 


3.72e-08 


5.81e-09 


2.16e-03 


1.40e+00 


N43 


1.03e-05 


1.64e-12 


9.89e-06 


2.65e-08 


4.14e-09 


1.33e-03 


1.69e+00 


N44 


7.00e-06 


1.05e-12 


4.74e-06 


1.38e-08 


2.27e-09 


8.56e-04 


1.47e+00 


N45 


8.27e-06 


1.36e-12 


5.99e-06 


1.58e-08 


2.62e-09 


8.19e-04 


1.73e+00 


N46 


6.15e-06 


1.21e-12 


4.08e-06 


l.lle-08 


1.85e-09 


6.41e-04 


1.68e+00 


N47 


5.77e-06 


1.02e-12 


2.85e-06 


l.lOe-08 


1.86e-09 


5.28e-04 


1.63e+00 


N48 


8.00e-06 


1.95e-12 


4.71e-06 


1.83e-08 


3.13e-09 


8.83e-04 


1.48e+00 


N49 


1.28e-05 


2.77e-12 


1.56e-05 


2.17e-08 


3.75e-09 


1.79e-03 


1.43e+00 


N50 


1.49e-05 


2.69e-12 


1.41e-05 


2.39e-08 


4.18e-09 


1.30e-03 


1.89e+00 


N51 


9.91e-06 


3.54e-12 


5.56e-06 


3.73e-08 


6.59e-09 


l.Ole-03 


1.96e+00 


N52 


6.98e-06 


3.22e-12 


2.96e-06 


3.51e-08 


6.27e-09 


8.13e-04 


1.83e+00 


N53 


4.48e-06 


1.31e-12 


1.37e-06 


1.41e-08 


2.56e-09 


4.78e-04 


1.42e+00 


N54 


3.72e-06 


1.12e-12 


2.21e-06 


6.24e-09 


1.14e-09 


3.22e-04 


1.78e+00 


N55 


8.24e-06 


1.76e-12 


1.43e-05 


1.34e-08 


2.48e-09 


8.23e-04 


1.72e+00 


N56 


5.23e-06 


1.96e-12 


1.26e-05 


1.50e-08 


2.80e-09 


7.96e-04 


1.69e+00 


N57-N60 


2.70e-06 


9.94e-13 


2.69e-06 


2.27e-09 


1.87e-10 


5.12e-04 


1.56e+00 


N61 


4.13e-06 


1.40e-12 


2.26e-06 


2.55e-09 


5.07e-10 


2.19e-04 


1.90e+00 


N62 


4.71e-06 


2.54e-12 


5.81e-06 


4.41e-09 


8.87e-10 


4.09e-04 


1.74e+00 


N63 


8.22e-06 


2.20e-12 


8.17e-06 


3.11e-09 


6.35e-10 


5.66e-04 


1.46e+00 


N65 


2.38e-06 


9.94e-13 


1.26e-06 


3.57e-10 


7.48e-ll 


1.47e-04 


1.75e+00 


N66 


4.10e-06 


8.42e-13 


8.15e-06 


3.86e-10 


9.97e-ll 


1.68e-04 


1.71e+00 


N67-N68 


2.72e-06 


1.34e-12 


4.44e-06 


4.26e-10 


3.10e-10 


1.89e-04 


2.70e+00 


N69 


2.76e-06 


1.20e-12 


1.31e-06 


4.89e-09 


3.45e-10 


8.09e-05 


2.08e+00 


N70-N72 


1.19e-06 


7.98e-13 


3.19e-07 


1.45e-09 


l.Ole-10 


5.06e-05 


2.12e+00 


N73-N77 


3.53e-06 


1.97e-12 


5.25e-07 


4.24e-10 


2.83e-ll 


8.51e-04 


2.31e+00 


given in erg s"' 


sr^'cm"^ at the angular resolution of the [C iil data. Values in 


bold face are upper limits 
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Fig. C.l: Examples of some of the [Cii] ISO/LWS spectra and SEDs obtained in this work. The vertical dashed lines mark the 
integration interval to derive the FIR continuum, 42.5 yum and 122.5jt/m ( Dale & Helou|2002 l. a,b: The first row shows the [Cii] 
spectrum and the SED of the inter-arm position S32. c,d: The second row shows the spectrum and SED from the nucleus of M 33. 
e,f: The third rows shows observations of the H ii region BCLMP 302 (N49). g,h: Finally, the last row shows the spectrum and tlkS 
SED from one of the positions (N61) in the outer, southern N2 region. The SEDs show MIPS 24/vm and 70jum, PACS 100/im and 
160/im and SPIRE 250/vm, 350jum and 500;um data. In the [C ii] spectra, green lines show Gaussian fits to the spectra and blue lines 
show fitted baselines. 
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Fig. D.l: [Cii] versus CO, normalized with the total infrared continuum (TIR). Big green filled circles show ISO/LWS data of the 
inner SI, Nl regions of M 33 while open circles show data of the outer S2, N2 regions. The lowest [C n]/CO ratio observed with 
ISO/LWS in M 33 is 1000 (lower blue dashed line), while the highest ratio is 41200 (upper blue dashed line). Black solid and dotted 
lines indicate lines of constant density n and FUV field Go, respectively, from the standard K99 PDR model with Ay - 10 mag and 
solar metallicity Z - I. 
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